Reevaluation of SNP heritability in complex human traits

SNP heritability, the proportion of phenotypic variance explained by SNPs, has been reported for many hundreds of traits. Its estimation requires strong prior assumptions about the distribution of heritability across the genome, but current assumptions have not been thoroughly tested. By analyzing imputed data for a large number of human traits, we empirically derive a model that more accurately describes how heritability varies with minor allele frequency (MAF), linkage disequilibrium (LD) and genotype certainty. Across 19 traits, our improved model leads to estimates of common SNP heritability on average 43% (s.d. 3%) higher than those obtained from the widely used software GCTA and 25% (s.d. 2%) higher than those from the recently proposed extension GCTA-LDMS. Previously, DNase I hypersensitivity sites were reported to explain 79% of SNP heritability; using our improved heritability model, their estimated contribution is only 24%

Triangulating evidence from observational and Mendelian randomization studies of ketone bodies for cognitive performance

BACKGROUND: Ketone bodies (KBs) are an alternative energy supply for brain functions when glucose is limited. The most abundant ketone metabolite, 3-β-hydroxybutyrate (BOHBUT), has been suggested to prevent or delay cognitive impairment, but the evidence remains unclear. We triangulated observational and Mendelian randomization (MR) studies to investigate the association and causation between KBs and cognitive function. METHODS: In observational analyses of 5506 participants aged ≥ 45 years from the Whitehall II study, we used multiple linear regression to investigate the associations between categorized KBs and cognitive function scores. Two-sample MR was carried out using summary statistics from an in-house KBs meta-analysis between the University College London-London School of Hygiene and Tropical Medicine-Edinburgh-Bristol (UCLEB) Consortium and Kettunen et al. (N = 45,031), and publicly available summary statistics of cognitive performance and Alzheimer's disease (AD) from the Social Science Genetic Association Consortium (N = 257,841), and the International Genomics of Alzheimer's Project (N = 54,162), respectively. Both strong (P < 5 × 10-8) and suggestive (P < 1 × 10-5) sets of instrumental variables for BOHBUT were applied. Finally, we performed cis-MR on OXCT1, a well-known gene for KB catabolism. RESULTS: BOHBUT was positively associated with general cognitive function (β = 0.26, P = 9.74 × 10-3). In MR analyses, we observed a protective effect of BOHBUT on cognitive performance (inverse variance weighted: βIVW = 7.89 × 10-2, PIVW = 1.03 × 10-2; weighted median: βW-Median = 8.65 × 10-2, PW-Median = 9.60 × 10-3) and a protective effect on AD (βIVW =  - 0.31, odds ratio: OR = 0.74, PIVW = 3.06 × 10-2). Cis-MR showed little evidence of therapeutic modulation of OXCT1 on cognitive impairment. CONCLUSIONS: Triangulation of evidence suggests that BOHBUT has a beneficial effect on cognitive performance. Our findings raise the hypothesis that increased BOHBUT may improve general cognitive functions, delaying cognitive impairment and reducing the risk of AD

Circulating Fatty Acids and Risk of Coronary Heart Disease and Stroke: Individual Participant Data Meta-Analysis in Up to 16 126 Participants

BACKGROUND We aimed at investigating the association of circulating fatty acids with coronary heart disease (CHD) and stroke risk. METHODS AND RESULTS We conducted an individual‐participant data meta‐analysis of 5 UK‐based cohorts and 1 matched case‐control study. Fatty acids (ie, omega‐3 docosahexaenoic acid, omega‐6 linoleic acid, monounsaturated and saturated fatty acids) were measured at baseline using an automated high‐throughput serum nuclear magnetic resonance metabolomics platform. Data from 3022 incident CHD cases (13 104 controls) and 1606 incident stroke cases (13 369 controls) were included. Logistic regression was used to model the relation between fatty acids and odds of CHD and stroke, adjusting for demographic and lifestyle variables only (ie, minimally adjusted model) or with further adjustment for other fatty acids (ie, fully adjusted model). Although circulating docosahexaenoic acid, but not linoleic acid, was related to lower CHD risk in the fully adjusted model (odds ratio, 0.85; 95% CI, 0.76–0.95 per standard unit of docosahexaenoic acid), there was evidence of high between‐study heterogeneity and effect modification by study design. Stroke risk was consistently lower with increasing circulating linoleic acid (odds ratio for fully adjusted model, 0.82; 95% CI, 0.75–0.90). Circulating monounsaturated fatty acids were associated with higher CHD risk across all models and with stroke risk in the fully adjusted model (odds ratio, 1.22; 95% CI, 1.03–1.44). Saturated fatty acids were not related to increased CHD risk in the fully adjusted model (odds ratio, 0.94; 95% CI, 0.82–1.09), or stroke risk. CONCLUSIONS We found consistent evidence that linoleic acid was associated with decreased risk of stroke and that monounsaturated fatty acids were associated with increased risk of CHD. The different pattern between CHD and stroke in terms of fatty acids risk profile suggests future studies should be cautious about using composite events. Different study designs are needed to assess which, if any, of the associations observed is causal

Variant rs10911021 that associates with coronary heart disease in type 2 diabetes, is associated with lower concentrations of circulating HDL cholesterol and large HDL particles but not with amino acids.

AIMS: An intergenic locus on chromosome 1 (lead SNP rs10911021) was previously associated with coronary heart disease (CHD) in type 2 diabetes (T2D). Using data from the UCLEB consortium we investigated the relationship between rs10911021 and CHD in T2D, whether rs10911021 was associated with levels of amino acids involved in the γ-glutamyl cycle or any conventional risk factors (CRFs) for CHD in the T2D participants. METHODS: Four UCLEB studies (n = 6531) had rs10911021 imputation, CHD in T2D, CRF and metabolomics data determined using a nuclear magnetic resonance based platform. RESULTS: The expected direction of effect between rs10911021 and CHD in T2D was observed (1377 no CHD/160 CHD; minor allele OR 0.80, 95 % CI 0.60-1.06) although this was not statistically significant (p = 0.13). No association between rs10911021 and CHD was seen in non-T2D participants (11218 no CHD/1274 CHD; minor allele OR 1.00 95 % CIs 0.92-1.10). In T2D participants, while no associations were observed between rs10911021 and the nine amino acids measured, rs10911021 was associated with HDL-cholesterol (p = 0.0005) but the minor "protective" allele was associated with lower levels (-0.034 mmol/l per allele). Focusing more closely on the HDL-cholesterol subclasses measured, we observed that rs10911021 was associated with six large HDL particle measures in T2D (all p < 0.001). No significant associations were seen in non-T2D subjects. CONCLUSIONS: Our findings are consistent with a true association between rs10911021 and CHD in T2D. The protective minor allele was associated with lower HDL-cholesterol and reductions in HDL particle traits. Our results indicate a complex relationship between rs10911021 and CHD in T2D

Establishing reference intervals for triglyceride containing lipoprotein sub-fraction metabolites measured using Nuclear Magnetic Resonance Spectroscopy in a UK population

Background Nuclear magnetic resonance (NMR) spectroscopy allows triglycerides to be subclassified into 14 different classes based on particle size and lipid content. We recently showed that these subfractions have differential associations with cardiovascular disease events. Here we report the distributions and define reference interval ranges for 14 triglyceride-containing lipoprotein subfraction metabolites. Methods Lipoprotein subfractions using the Nightingale NMR platform were measured in 9073 participants from four cohort studies contributing to the UCL-Edinburgh-Bristol consortium. The distribution of each metabolite was assessed, and reference interval ranges were calculated for a disease-free population, by sex and age group (65 years), and in a subgroup population of participants with cardiovascular disease or type 2 diabetes. We also determined the distribution across body mass index and smoking status. Results The largest reference interval range was observed in the medium very-low density lipoprotein subclass (2.5th 97.5th percentile; 0.08 to 0.68 mmol/L). The reference intervals were comparable among male and female participants, with the exception of triglyceride in high-density lipoprotein. Triglyceride subfraction concentrations in very-low density lipoprotein, intermediate-density lipoprotein, low-density lipoprotein and high-density lipoprotein subclasses increased with increasing age and increasing body mass index. Triglyceride subfraction concentrations were significantly higher in ever smokers compared to never smokers, among those with clinical chemistry measured total triglyceride greater than 1.7 mmol/L, and in those with cardiovascular disease, and type 2 diabetes as compared to disease-free subjects. Conclusion This is the first study to establish reference interval ranges for 14 triglyceride-containing lipoprotein subfractions in samples from the general population measured using the nuclear magnetic resonance platform. The utility of nuclear magnetic resonance lipid measures may lead to greater insights for the role of triglyceride in cardiovascular disease, emphasizing the importance of appropriate reference interval ranges for future clinical decision making

Sixty-five common genetic variants and prediction of type 2 diabetes.

We developed a 65 type 2 diabetes (T2D) variant-weighted gene score to examine the impact on T2D risk assessment in a U.K.-based consortium of prospective studies, with subjects initially free from T2D (N = 13,294; 37.3% women; mean age 58.5 [38-99] years). We compared the performance of the gene score with the phenotypically derived Framingham Offspring Study T2D risk model and then the two in combination. Over the median 10 years of follow-up, 804 participants developed T2D. The odds ratio for T2D (top vs. bottom quintiles of gene score) was 2.70 (95% CI 2.12-3.43). With a 10% false-positive rate, the genetic score alone detected 19.9% incident cases, the Framingham risk model 30.7%, and together 37.3%. The respective area under the receiver operator characteristic curves were 0.60 (95% CI 0.58-0.62), 0.75 (95% CI 0.73 to 0.77), and 0.76 (95% CI 0.75 to 0.78). The combined risk score net reclassification improvement (NRI) was 8.1% (5.0 to 11.2; P = 3.31 × 10(-7)). While BMI stratification into tertiles influenced the NRI (BMI ≤24.5 kg/m(2), 27.6% [95% CI 17.7-37.5], P = 4.82 × 10(-8); 24.5-27.5 kg/m(2), 11.6% [95% CI 5.8-17.4], P = 9.88 × 10(-5); >27.5 kg/m(2), 2.6% [95% CI -1.4 to 6.6], P = 0.20), age categories did not. The addition of the gene score to a phenotypic risk model leads to a potentially clinically important improvement in discrimination of incident T2D

Functional Analysis of the Coronary Heart Disease Risk Locus on Chromosome 21q22.

Background. The coronary heart disease (CHD) risk locus on 21q22 (lead SNP rs9982601) lies within a "gene desert." The aim of this study was to assess if this locus is associated with CHD risk factors and to identify the functional variant(s) and gene(s) involved. Methods. A phenome scan was performed with UCLEB Consortium data. Allele-specific protein binding was studied using electrophoretic mobility shift assays. Dual-reporter luciferase assays were used to assess the impact of genetic variation on expression. Expression quantitative trait analysis was performed with Advanced Study of Aortic Pathology (ASAP) and Genotype-Tissue Expression (GTEx) consortium data. Results. A suggestive association between QT interval and the locus was observed (rs9982601  p = 0.04). One variant at the locus, rs28451064, showed allele-specific protein binding and its minor allele showed 12% higher luciferase expression (p = 4.82 × 10-3) compared to the common allele. The minor allele of rs9982601 was associated with higher expression of the closest upstream genes (SLC5A3 1.30-fold increase p = 3.98 × 10-5; MRPS6 1.15-fold increase p = 9.60 × 10-4) in aortic intima media in ASAP. Both rs9982601 and rs28451064 showed a suggestive association with MRPS6 expression in relevant tissues in the GTEx data. Conclusions. A candidate functional variant, rs28451064, was identified. Future work should focus on identifying the pathway(s) involved

Plasma urate concentration and risk of coronary heart disease: a Mendelian randomisation analysis.

BACKGROUND: Increased circulating plasma urate concentration is associated with an increased risk of coronary heart disease, but the extent of any causative effect of urate on risk of coronary heart disease is still unclear. In this study, we aimed to clarify any causal role of urate on coronary heart disease risk using Mendelian randomisation analysis. METHODS: We first did a fixed-effects meta-analysis of the observational association of plasma urate and risk of coronary heart disease. We then used a conventional Mendelian randomisation approach to investigate the causal relevance using a genetic instrument based on 31 urate-associated single nucleotide polymorphisms (SNPs). To account for potential pleiotropic associations of certain SNPs with risk factors other than urate, we additionally did both a multivariable Mendelian randomisation analysis, in which the genetic associations of SNPs with systolic and diastolic blood pressure, HDL cholesterol, and triglycerides were included as covariates, and an Egger Mendelian randomisation (MR-Egger) analysis to estimate a causal effect accounting for unmeasured pleiotropy. FINDINGS: In the meta-analysis of 17 prospective observational studies (166 486 individuals; 9784 coronary heart disease events) a 1 SD higher urate concentration was associated with an odds ratio (OR) for coronary heart disease of 1·07 (95% CI 1·04-1·10). The corresponding OR estimates from the conventional, multivariable adjusted, and Egger Mendelian randomisation analysis (58 studies; 198 598 individuals; 65 877 events) were 1·18 (95% CI 1·08-1·29), 1·10 (1·00-1·22), and 1·05 (0·92-1·20), respectively, per 1 SD increment in plasma urate. INTERPRETATION: Conventional and multivariate Mendelian randomisation analysis implicates a causal role for urate in the development of coronary heart disease, but these estimates might be inflated by hidden pleiotropy. Egger Mendelian randomisation analysis, which accounts for pleiotropy but has less statistical power, suggests there might be no causal effect. These results might help investigators to determine the priority of trials of urate lowering for the prevention of coronary heart disease compared with other potential interventions. FUNDING: UK National Institute for Health Research, British Heart Foundation, and UK Medical Research Council

PCSK9 genetic variants and risk of type 2 diabetes: a mendelian randomisation study

BACKGROUND: Statin treatment and variants in the gene encoding HMG-CoA reductase are associated with reductions in both the concentration of LDL cholesterol and the risk of coronary heart disease, but also with modest hyperglycaemia, increased bodyweight, and modestly increased risk of type 2 diabetes, which in no way offsets their substantial benefits. We sought to investigate the associations of LDL cholesterol-lowering PCSK9 variants with type 2 diabetes and related biomarkers to gauge the likely effects of PCSK9 inhibitors on diabetes risk. METHODS: In this mendelian randomisation study, we used data from cohort studies, randomised controlled trials, case control studies, and genetic consortia to estimate associations of PCSK9 genetic variants with LDL cholesterol, fasting blood glucose, HbA1c, fasting insulin, bodyweight, waist-to-hip ratio, BMI, and risk of type 2 diabetes, using a standardised analysis plan, meta-analyses, and weighted gene-centric scores. FINDINGS: Data were available for more than 550 000 individuals and 51 623 cases of type 2 diabetes. Combined analyses of four independent PCSK9 variants (rs11583680, rs11591147, rs2479409, and rs11206510) scaled to 1 mmol/L lower LDL cholesterol showed associations with increased fasting glucose (0·09 mmol/L, 95% CI 0·02 to 0·15), bodyweight (1·03 kg, 0·24 to 1·82), waist-to-hip ratio (0·006, 0·003 to 0·010), and an odds ratio for type diabetes of 1·29 (1·11 to 1·50). Based on the collected data, we did not identify associations with HbA1c (0·03%, -0·01 to 0·08), fasting insulin (0·00%, -0·06 to 0·07), and BMI (0·11 kg/m(2), -0·09 to 0·30). INTERPRETATION: PCSK9 variants associated with lower LDL cholesterol were also associated with circulating higher fasting glucose concentration, bodyweight, and waist-to-hip ratio, and an increased risk of type 2 diabetes. In trials of PCSK9 inhibitor drugs, investigators should carefully assess these safety outcomes and quantify the risks and benefits of PCSK9 inhibitor treatment, as was previously done for statins. FUNDING: British Heart Foundation, and University College London Hospitals NHS Foundation Trust (UCLH) National Institute for Health Research (NIHR) Biomedical Research Centre.This work was supported by a British Heart Foundation Programme Grant (RG/10/12/28456). AFS is funded by University College London Hospitals NHS Foundation Trust (UCLH) National Institute for Health Research (NIHR) Biomedical Research Centre (BRC10200) and by a UCL springboard population science fellowship. FWA is supported by a Dekker scholarship-Junior Staff Member 2014T001–Netherlands Heart Foundation and UCL Hospitals NIHR Biomedical Research Centre. ADH is an NIHR Senior Investigator. Funding information and acknowledgments for studies contributing data are reported in the appendix

The UK10K project identifies rare variants in health and disease

M. Kivimäki työryhmäjäsen.The contribution of rare and low-frequency variants to human traits is largely unexplored. Here we describe insights from sequencing whole genomes (low read depth, 7x) or exomes (high read depth, 80x) of nearly 10,000 individuals from population-based and disease collections. In extensively phenotyped cohorts we characterize over 24 million novel sequence variants, generate a highly accurate imputation reference panel and identify novel alleles associated with levels of triglycerides (APOB), adiponectin (ADIPOQ) and low-density lipoprotein cholesterol (LDLR and RGAG1) from single-marker and rare variant aggregation tests. We describe population structure and functional annotation of rare and low-frequency variants, use the data to estimate the benefits of sequencing for association studies, and summarize lessons from disease-specific collections. Finally, we make available an extensive resource, including individual-level genetic and phenotypic data and web-based tools to facilitate the exploration of association results.Peer reviewe